Statistical theory of quasi stationary states beyond the single water-bag case study 
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An analytical solution for the out-of-equilibrium quasi-stationary states of the paradigmatic 
Hamiltonian Mean Field (HMF) model can be obtained from a maximum entropy principle. The 
theory has been so far tested with reference to a specific class of initial condition, the so called 
(single-level) water-bag type. In this paper a step forward is taken by considering an arbitrary 
number of overlapping water bags. The theory is benchmarked to direct microcanonical simulations 
performed for the case of a two-levels water-bag. The comparison is shown to return an excellent 
agreement. 



PACS numbers: 05.20.-y, 05.45.-a, 05.70.Ce, 05.70.Fh 

I. INTRODUCTION 

Long range interacting systems (LRS) are becoming a 
popular topic of investigation Q in physics, due to the 
rich and intriguing phenomenology that they display. A 
system is said to fall in the realm of LRS if the two body 
potential scales as r~ a with a < d, and where r stands 
for the inter-particle distance and d the dimension of the 
embedding space. Several physical systems share this 
property, which ideally embraces distinct domains of ap- 
plications. Gravity [3| is certainly the most spectacular 
example among the wide gallery of systems governed by 
long range interactions, but equally important are the 
cases of turbulence Q , plasmas Q and wave- particle in- 
teractions 

Peculiar and counterintuitive thermodynamics features 
manifest in LRS: negative specific heat can occasionally 
develop in the microcanonical ensemble, close to first or- 
der phase transitions, a surprising fact first discovered in 
astrophysical context, that seeds statistical ensemble in- 
equivalence. As concerns the dynamics, LRS have been 
reported to experience a very slow relaxation towards 
the deputed thermodynamic equilibrium. Indeed, they 
can be trapped in long lasting out-of-equilibrium phases 
called Quasi Stationary States (QSSs). The lifetime of 
the QSSs diverges with the systems size N. Interestingly, 
it displays different scaling behaviors versus N, which 
range from exponential to power law, being relic of the 
specific initial condition selected. As a consequence, the 
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orders the limits N — > oo and t — > oo are taken do mat- 
ter. Performing the continuum limit before the infinite 
time limit, implies preventing the system from eventually 
attaining its equilibrium and so freezing it indefinitely in 
the QSS phase. 

In physical applications where long-range couplings are 
at play, as the ones mentioned above, the number of ele- 
mentary constituents composing the system being exam- 
ined, is generally large. The time of duration of the out- 
of-equilibrium phase can therefore be exceedingly long, 
definitely longer the the time of observation to which ex- 
perimentalists are bound. Given this scenario, it is of 
paramount importance to develop dedicated analytical 
strategies to gain quantitative insight into the complex 
and diverse zoology of the QSSs, as revealed by direct 
numerical simulations. Working along these lines, it was 
shown that QSS can be successfully interpreted as equi- 
libria of the collisionless Vlasov equation which appears 
to rule the dynamics of a broad family of long range mod- 
els, when recovering the continuum picture from the gov- 
erning discrete formulation. The average characteristics 
of the QSS, including the emergence of out-of-equilibrium 
transitions, can be analytically predicted via a maxi- 
mum entropy variational principle, pioneered by Lynden- 
Bell in [l4[ and more recently revisited with reference to 
paradigmatic long-range applications @, H, HU ■ 
As we shall clarify in the forthcoming sections, the pre- 
dictive adequacy of the Lyndcn-Bcll violent relaxation 
theory has been so far solely assessed for a very specific 
class of initial conditions. These are the so called (single) 
water-bags: particles are assumed to initially populate 
a bound domain of phase space and therein distributed 
with a uniform probability. The aim of this paper is to 
take one simple step forward and challenge the validity of 
the theory when particles are instead distributed within 
two (uniformly filled) levels. In principle, any smooth 
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profile could be approximated by a pieccwise function, 
made of an arbitrary number of collated water- bags fl6| . 
Our idea is to perform a first step towards the general- 
ized multi-levels setting, by first evaluating the formal 
complexity of the procedure involved and then drawing a 
direct comparison with the numerics relative to the two 
level case. 

To accomplish this task, we will focus on the cele- 
brated Hamiltonian Mean Field model, often referred to 
as the representative model of long range interactions. 
The HMF describes the motion on a circle of an ensem- 
ble of N rotors mutually coupled via an all-to-all cosines 
like potential. In the continuum limit, the single parti- 
cle distribution function obeys to the Vlasov equation, 
the driving potential being self-consistently provided by 
the global magnetization, namely the degree of inher- 
ent bunching. QSSs exist for the HMF model and have 
been deeply studied, both with analytical and numerical 
means, for the single water-bag case. 

The paper is organized as follows. In the next section 
we shall introduce the discrete HMF model and discuss 
its continuum Vlasov based representation. We will also 
introduce the basic of the violent relaxation approach. 
Then, in section III, we shall turn to discussing the gen- 
eralized water bag setting, solving, in section IV, the cor- 
responding variational problem. We will then specialize 
in section V on the two levels case and compare the the- 
ory predictions to the simulations. Finally in section VI 
we will sum up and conclude. 

II. THE HMF MODEL 

The HMF model describes the dynamics of N particles 
(rotors) moving on a circle and interacting via a mean 
field potential which is self-consistently generated by the 
particles themselves. Formally, the HMF is defined by 
the following Hamiltonian: 

N 2 i N 

i—l ij — 1 

where (9; identifies the position of particle i on the 
circle and pi is the canonically conjugated momentum. 
Because the interactions are not bounded to a small 
number of neighboring particles, the interaction is 
all-to-all thus the potential is inherently long range. 

Starting from the water bag initial condition, the HMF 
system experiences a fast relaxation towards an interme- 
diate regime, before the final equilibrium is eventually 
attained. This metastable phase is a Quasi Stationary 
State (QSS), the out-of-equilibrium transient to which 
we have alluded to in the introduction. The lifetime of 
the QSS is shown to diverge with the system size N, an 
observation that has non-trivial consequences when one 
wishes to inspect the continuum (N oo) limit. QSSs 
are in fact stable, attractive equilibria of the continuous 



analogue of the discrete Hamiltonian picture, and bear 
distinctive traits that make them substantially different 
from the corresponding equilibrium solutions. 

To monitor the dynamics of the system it is customary 
to record the time evolution of the magnetization. This 
latter is defined as: 

N 

M = ^y e . (2) 

i=l 

It is a complex quantity whose modulus M measures 
the degree of bunching of the distribution of particles. 
Depending on the selected characteristics of the initial 
(single) water-bag, the system can evolve towards an (al- 
most) homogeneous QSS or, conversely, result in a mag- 
netized phase. The swap between the two regimes can 
be understood as a genuine phase transition, with en- 
ergy and initial magnetization playing the role of control 
parameters. 

It can be rigorously shown [l3j that, in the contin- 
uum limit, the HMF system is formally described by the 
Vlasov equation, which governs the evolution of the sin- 
gle particle distribution function f{9,p 1 t). In formula: 

%+P%- ( M * [/] sin 6 ~ M v I/I cos °) |£ = °. ( 3 ) 

where M x — J f cos QdQdp and M y = J f sin OdOdp are 
the two components of the magnetization M. With refer- 
ence to cosmological applications, Lyndcn-Bcll proposed 
an analytical approach to determine the stationary solu- 
tions of the Vlasov equation, pioneering the theory that 
it is nowadays referred to as to the violent relaxation the- 
ory. He first considered the coarse grained distribution 
/, obtained by averaging the microscopic f(6,p,t) over 
a finite grid. Then, the key idea is to associate to / a 
mixing entropy S[f], via a rigorous counting of the mi- 
croscopic configurations that are compatible with a given 
macroscopic state. The steps involved in the derivation 
are highlighted in the remaining part of this section. 

Let us label with p(9,p,rj) the probability density of 
finding the level of phase density 77 in the neighborhood of 
the position [0,p] in phase space. We are here implicitly 
assuming to deal with a continuum spectrum of allowed 
levels. Following Lynden-Bell the coarse grained locally 
averaged single particle distribution reads: 

f(6,p)= J p(6,p,ri)ridr,. (4) 

The Vlasov equation, which rules the dynamics of /, 
conserves the hypervolumes 1/(77) associated to each of 
the selected levels [lT}. Mathematically, the quantities 
1/(77) reads: 

v{rj) = [ p(e,p,r])dedp , (5) 



and arc therefore invariant of the dynamics. Assume now 
to deal with a discrete set of n distinct levels. The prob- 
ability density function reads therefore: 



Pn (0,p,v) = ^2f(9,P) a i$(v - fi) + 

i=l 

(l-£/(0,p)a^ <%), 



(6) 



where fi refers to the value of the «-th level and a,; stands 
for the portion of phase space that is hosting the se- 
lected level. The rightmost term in the previous equa- 
tion, stands thus for the background. We shall return on 
these aspects later on, when commenting on the under- 
lying normalization condition. 

Starting from this setting, Lynden-Bell suggested to 
divide the phase space macrocell [9, 9 + d9;p,p + dp] into 
microcells each occupied by just one of the allowed lev- 
els. The Lynden-Bell mixing entropy is obtained as the 
logarithm of the total number of microstates associated 
to a given macrostatc. A combinatorial calculation [l?} 
yields to: 



S 



Pn (9, p, t]) In p n (9, p, rj) dOdpdrj, 



(7) 



The statistical equilibrium of the system follows from 
maximizing the entropy functional S while imposing the 
constraints of the dynamics: energy, momentum and nor- 
malization arc in fact conserved quantities, as well as the 
quantities v(rj). In the following we shall discuss a spe- 
cific class of initial condition, the multi-level watcrbags, 
which naturally extends beyond the single water bag case 
study, so far explicitly considered in the literature. It is 
our intention to test the predictive ability of the Lynden- 
Bell theory within such generalized framework. The the- 
ory will be developed with reference to the general set- 
ting, including n levels. The benchmark with direct sim- 
ulations will be instead limited to the two-levels case, i.e. 
n = 2. 




FIG. 1: Pictorial representation of a tree-levels (n = 3) water- 
bag initial condition. 



a Fcrmionic stationary distribution, which successfully 
enables to capture some of the essential traits of the 
QSS. These includes an accurate characterization of the 
out-of-equilibrium transitions from magnetized to non- 
magnetized QSS. First and second order phase transi- 
tions, that merge in a tricritical point, were in fact sin- 
gled out for the HMF model, a theoretical prediction con- 
firmed by direct numerical inspection. As stated above, 
the general philosophy that inspires the Lynden-Bell the- 
ory is however broader than the specific realm to which it 
was relegated and its potentiality deserves to be further 
clarified. We will here extend the treatment to the multi- 
levels water bag initial condition, a step that opens up 
the perspective to eventually handle more realistic sce- 
narios, where smooth distributions could be considered. 

Following the notation introduced above, the arbitrary 
integer n quantifies the total number of distinct levels 
that are to be allowed for, when considering the general- 
ized initial distribution function fi n u- Arguably, by ac- 
counting for a large enough collection of independent and 
discrete levels, one can approximately mimic any smooth 
profile. A pictorial representation of the family of initial 
conditions to which we shall refer to in the forthcoming 
sections when discussing the specific case study n = 3 is 
depicted in figure Q] 

Mathematically, the initial distribution function fi n u 
can be written as: 



,P) = 



fi 



if 6 S 0j and p G Pj, 



/o = elsewhere, 



(8) 



III. THE GENERALIZED WATER-BAG 

The single water-bag initial condition takes a constant 
value /i within a finite portion of the phase space, and 
zero outside of it. Although this is the only prescrip- 
tion to be accommodated for, rectangular domains are 
usually chosen for practical computational reasons. Fol- 
lowing [l5j], we shall label [A0, Ap] the widths of such 
a rectangle, as calculated respectively along 9 and p di- 
rections. A second simplification is also customarily be 
assumed: the rectangle is centered in the origin, so that 

By operating in this context, the Lynden-Bell varia- 
tional problem studied in e.g. [l5| is shown to yield to 



Here Tj = [Qj,Pj], j — l,...,n identifies the domain 
in phase space associated to level fj. The corresponding 
volume is labeled ctj. 

A simple algebraic manipulation starting from the def- 
inition of the probability density ([5]) yields to: 

/ finit(9,p)d9dp = / p n (9,p,r])r]d9dpdr] = ^fjatj = 1. 

J J j=0 

(9) 

The scalar relation ^ links together the 2n constants, 
fj and dj, that are to be assigned to fully specify the 
initial condition. In other words, only 2n — 1 scalars are 
needed to completely parameterize the initial condition. 
Importantly, the single water bag limit is readily recov- 
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ered once the phase space support of the levels indexed 
with j other than j ' = 1 shrinks and eventually fades out. 
This condition implies requiring otj — > for j > 1. In 
formula: 



Pi 



lim p n 

->0,j>l 



f ai 6( v -A) + (l-/«i) 5( V ) 



(10) 

Moreover, by making use of the normalization condition 
one gets ct\ = 4-, which, inserted (|10[) . returns immedi- 
ately the well known form of the one level density dis- 
tribution function. Notice that the initial value of the 
macroscopic observables can be computed from the ex- 
plicit knowledge of p\ . 



IV. THE GENERALIZED N-LEVELS 
EQUILIBRIUM 



j = 1, ..n, defined as: 



f(9,p)a j d6dp, 



(14) 



Moreover using the normalization condition for the 
coarse grained distribution function f(9,p), we get 
fj [/] = aj . Equivalently, by imposing the above con- 
straints on the hypcrvolumcs, we also guarantee the nor- 
malization of the distribution function, which physically 
amounts to impose the conservation of the mass. 

Summing up, the variational problem that needs to be 
solved to eventually recover the stationary distribution 
f eq {0,p) reads: 

max{5[/] \E[f]= E n ;P [f] = P n ; Vl [f] = a,}, (15) 

where the entropy functional S[f] is given by eq. (|11[) . 
This immediately translates into: 



The QSS distribution function f eq (9,p) for the HMF 
model, relative to the generalized n-levels water-bag 
initial condition, is found by maximizing the Lynden 
Bell entropy, under the constrains of the dynamics. 
This in turn implies solving a variational problem. The 
solution is relative to the microcanonical ensemble, since 
the Vlasov equation implies that we work with fixed 
total energy. 

From equation ([7]) and ©, the generic n-levels entropy 
takes the following functional form: 



/n 

3 = 1 



(11) 



(i - y, ln (i - E d9d p 

3=1 3=1 



The conserved quantities are respectively the energy E: 

M[jf - 1 



E U\= J \f{0,P)d8d V - 



and the total momentum P 



P[f] 



f(6,p)pd6dp = P n 



E n , (12) 



(13) 



The scalar quantity E n relates to the geometric char- 
acteristics of the bounded domains that define our 
initial condition. Conversely, as we will be dealing with 
patches Tj symmetric with respect to the origin, one can 
immediately realize that P n = 0. 

The n volumes of phase space, each deputed to host- 
ing one of the considered levels, are also invariant of the 
dynamics. We have therefore to account for the conser- 
vation of n additional quantities, the volumes Vj [/] for 



SS - I35E - XSP - Y N Sv 3 = °- (16) 

3 = 1 

where /3, A and stands for the Lagrange multipliers 
associated respectively to energy, momentum and vol- 
umes (or equivalently mass) conservations. 

A straightforward calculation yields to the following 
expression for f eq (d,p)' 



feq 



1 



/3'(4-M[/ c ,].m)+A' 



where 



where 



B + Ae 



B = Y a >>A= n«? 

3=1 \3=1 



P B ' 
A'- A 



(17) 



(18) 



(19) 



En 
3=1 »i 



B 



and m = [cos(#), sin(#)]. 

The above solution is clearly consistent with that ob- 
tained for the single water bag case study . This latter 
is in fact recovered in the limit {aj — > for j > 1 while 
at = i): 



lim f eq = 

aj ^o,j>i 



fi 



(20) 
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Notice that the equilibrium distribution f eq depends on 
M, which is in turn function of f eq itself. The two com- 
ponents of the magnetization, respectively M x and M y 
are therefore unknowns of the problem, implicitly depen- 
dent on f eq . This latter is parameterized in terms of 
the Lagrange multipliers. Their values need to be self- 
consistcntly singled out. As a first simplification, we ob- 
serve that the specific symmetry of the selected initial 
condition (P n = 0) implies A = 0. Hence, just the two 
residual Lagrange multipliers are to be computed: the 
Lynden-Bcll inverse temperature (3 and the cumulative 
chemical potential // fl9j . The number of unknowns total 
therefore to four (M x , M y , (3, //) and enter the following 
system of implicit equations for the constraints: 



E = — i- / e^ m F 2 (y) d9 - 

2/3' 3 / 2 J yyi 2 

l = -^je^™F {y)d6 

M x = ^= J eP' M m F Q (y)cos(e)d9 

M y = -^=J e^' M m F (y) sin(0) d0 



(21) 
(22) 
(23) 
(24) 



Here we have expressed the relations as function 

of the Fermi integrals Fh (y) = / f + e *Ji /2 dp, with 

y = ABe^' m ' m and A = A^ x e~^' . The system of equa- 
tions pljl .([22 ]l .p3 ]) . (|2l]l can be solved numerically. In 
doing so one obtains a numerical value for the involved 
Lagrange multipliers, as well as for the magnetization 
components, by varying the parameters that encode for 
the initial condition. We numerically checked (data not 
shown) that in the limit of a single water bag tXj>i — > 
the solution reported in [ll| is indeed recovered. In 
the following section we turn to discussing the theory 
predictions with reference to the simple case of two 
water bag (n = 2), validating the results versus direct 
numerical simulations. 



V. THE CASE n = 2: THEORY PREDICTIONS 
AND NUMERICAL SIMULATIONS. 

We here consider the simplifying setting where two 
levels (n = 2) water-bag are allowed for. We are in 
particular interested in monitoring the dependence of 

M = J~M% + versus the various parameters that 

characterize the initial condition. We recall in fact that, 
for the case of a single water-bag, out of equilibrium tran- 
sitions have been found 15| , which separates between ho- 
mogeneous and magnetized phases. A natural question 
is thus to understand what is going to happen if one ad- 



ditional level is introduced in the initial condition. The 
level fi is associated to a rectangular domain Ti of re- 
spective widths A6i and Api . The level / 2 insists instead 
on an adjacent domain r 2 , whose external perimeter is 
delimited by a rectangle of dimensions A9 2 and Ap 2 ■ The 
corresponding surface totals hence AQ 2 Ap 2 — A#iApi. 

Recall that the energy E2 (E n for n = 2) can be es- 
timated as dictated by formula (Ti"2"|) and reads in this 
specific case: 



E 2 = -(f 1 A6 1 Ap 3 1 + (f 2 -f l )f 2 A9 2 Apl) 



(25) 



1 - 16(/iApi sinA0i/2 + {h - /i)Apa sinA0 2 /2) 



The one-level limit is readily recovered by simultane- 
ously imposing A6* 2 — > and Ap 2 — > (which also im- 
plies a 2 — > 0). By invoking the normalization condition 
© the following relation holds: 



lim 

A0 2 ,Ap 2 



E 2 = \Ap\+ l -{l-Ml) (26) 
>0 o 2 



where Mq = 2sin(A0i/2)/A#i. The above relation coin- 
cides with the canonical expression for E\ , as e.g. derived 
in 0. 

Relation ([25]) enables us to estimate the energy asso- 
ciated to the selected initial condition and can be used 
in the self-consistency equations (|21[) . Before turning to 
illustrate the predicted solution, let us note that the nor- 
malization © reduces for n = 2 to: 



+ "2/2 = !• 



(27) 



To explore the parameter space we have decided to 
monitor the dependence of Mon/i, which therefore acts 
as a control parameter. To this end, we proceed by fixing 
the quantity A/ = / 2 — /1, the difference in hight of the 
considered levels. Furthermore, we specify the quantity 
ot\, while oi2 is calculated so to match the normalization 
constraint. 

The analysis is then repeated for distinct choices of 
A/, so to eventually elaborate on the importance of such 
crucial parameter. The results are displayed in figure [5] 
The curves collapse towards a point that corresponds to 
the limiting condition a 2 — > (fi = 1/ai): this spe- 
cial solution is met when the hypervolume populated by 
the level / 2 shrinks to zero, so driving the system to- 
wards the standard one level setting. By progressively 
reducing f\ the predicted magnetization first increases 
and subsequently decreases to eventually reach zero at a 
critical threshold /j\ For fi > /f the system is predicted 
to evolve towards a magnetized, hence non homogeneous 
phase. Alternatively, for < /1 < /[ a homogeneous 
phase is expected to occur. Interestingly, the transition 
point fi depends on the selected A/: the larger Af the 
smaller the value of the transition point, corresponding 



6 




Fast relaxation 
QSS 



Thermal equilibrium 
state 




FIG. 2: Analytical predictions for the equilibrium magnetiza- 
tion M[f eq ] as obtained for different values of the initial two 
levels water-bag distribution. The two levels are respectively 
labeled fx and /a. We here work at constant ax = 5 and 
Af = /a — /i, while moving the control parameter /i. The 
analysis is repeated for distinct values of Af (from left to right 
Af = 0.2,0.15,0.1,0.05 ). a.2 is computed according to eq. 
©. For f\ — ¥ 1/ax = 0.2 the normalization condition yields 
to Q2 — ¥ 0, and the distribution collapses to the limiting case 
of a single water-bag. 




FIG. 3: The analytical curves (same setting as in figure ([2J) 
are now plotted in the plane (E2, fi). (from left to right Af = 
0.2,0.15,0.1,0.05 ). We here only represent the points that 
are associated to positive M[f eq ]. The transition occurs at 
constant energy E2 — 0.675, regardless of the specific domains 
that result in the two level water-bag distribution. 



to a shift to the left in figure [5J Notice that above a 
limiting value of A/, which self-consistently corresponds 
to imposing a 2 > A/, the value of fx has to forcefully 
become negative so to respect the normalization condi- 
tion. A well hence opens up in phase space, an intriguing 
scenario that can be formally handled within the descrip- 
tive Vlasov framework but that we have here deliberately 
omitted to deepen any further. The smooth phase tran- 
sition as depicted in figure [5] is therefore lost above a 
threshold value of A/, when the predicted value of M 
associated to fx = turns out to be greater than zero. 

To elucidate the specificity of the outlined transition, 
we plot in figure[3]the energy E 2 , associated to each of the 
selected initial conditions, versus fx, for the same selec- 
tion of parameters as employed in figure [2j As suggested 
by visual inspection of the figure, the transitions, which 
we recall take place within a finite window in fx, always 
occur for an identical value of the energy (in this case 



FIG. 4: Magnetization M as a function of time t, as seen in a 
typical simulation. The system experiences a fast growth and 
then settle down into the lethargic QSS phase, whose dura- 
tion (data not shown) increases with N. Later on the system 
moves towards the deputed equilibrium. In this simulation 
a two levels water-bag is assumed with fx = 0.14, fi — 0.1, 
ax = 0.2, q 2 = 0.3. The energy is E 2 = 1.0 and N = 10 4 . 



E% — 0.675). The transition point is hence insensitive to 
the specificity of the two water-bags, being neither de- 
pendent on their associated volumes nor relative heights. 
It is in principle possible to extend the above analysis 
and so reconstruct the complete transition surface in the 
(/i) ii-, E2) space, a task which proves however demand- 
ing from the computational viewpoint and falls outside 
the scope of the present paper. 

To test the validity of the theory we have run a se- 
ries of numerical simulations of the HMF model. The 
implementation is based on fifth order McLachlan-Atela 
algorithm [l8[ with a time-step St — 0.1. The initial con- 
dition is of a two levels water-bag type, with respective 
domains assigned as follows the aforementioned prescrip- 
tions. As a preliminary check we have monitored the 
approach to equilibrium, figure 2] 

As expected and generalizing the conclusion that have 
been shown to hold for the simpler one level water-bag 
family of initial conditions, the system settles down into 
a QSS, whose lifetime grows with the number of simu- 
lated particles (data not shown). The QSS are indeed 
the target of our analysis and it is the magnetization 
as recorded in the QSS phase that needs to be compared 
to the Lynden-Bcll predictions. The comparison between 
theory and simulations is reported in figure[5J Filled sym- 
bols refer to the simulation while the solid line stand for 
the theory, for two distinct choice of Af. The agreement 
is certainly satisfying and points to the validity of the 
Lynden-Bell interpretative framework, beyond the case 
of the single water-bag, so far discussed in the literature. 



VI. CONCLUSIONS 

The dynamics of long range interacting system is stud- 
ied, as concerns the intriguing emergence of long lasting 
Quasi Stationary States. The problem is tackled within 
the context of the Hamiltonian Mean Field model, a very 
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FIG. 5: The analytical predictions (solid line) for the QSS 
magnetization as a function of /i in the two-levels water-bag 
case, are compared (filled circles) to the numerical simulations 
performed for N = 10 4 . The comparison is drawn for two 
distinct values of A/ (A/ = 0.1 (left) and A/ = 0.2 (right)). 
oti = 5 and 0:2 follows the normalization condition. Numerical 
values of M are computed as a time average over a finite time 
window where the QSS holds. The data are further mediated 
over 4 independent realizations. Expected uncertainties are 
about the size of the circle. 



popular and paradigmatic case study. Building on previ- 
ous evidences, the QSS are interpreted as stable equilib- 
ria of the Vlasov equation, which rules the dynamics of 
the discrete HMF system in the infinite system size limit 
(N — > 00). The QSS are hence characterized analytically 
by means of a maximum entropy principle inspired to the 
seminal work of Lyndcn Bell. This technique is known 
to yield to reliable predictions, when dealing with a very 
specific class of initial condition, the so called (single) 
water-bag. The scope of this paper is to push forward 
the analysis by considering the case where multiple wa- 
ter bags are allowed for. The theory is challenged with 
reference to the case of a two levels water-bag initial con- 
dition and the comparison with the simulations proves 
accurate. Phase transitions are in fact predicted and ob- 
served in direct ./V-body simulations. Motivated by this 
success, we argue that the Lyndcn-Bcll approach could 
be adapted to more complex, and so realistic, family of 
initial conditions. 
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